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RIGIDITY AND SHARP STABILITY ESTIMATES FOR 
HYPERSURFACES WITH CONSTANT AND ALMOST-CONSTANT 
NONLOCAL MEAN CURVATURE 

G. CIRAOLO, A. FIGALLI, F. MAGGI, AND M. NOVAGA 


Abstract. We prove that the boundary of a (not necessarily connected) bounded smooth set 
with constant nonlocal mean curvature is a sphere. More generally, and in contrast with what 
happens in the classical case, we show that the Lipschitz constant of the nonlocal mean curvature 
of such a boundary controls its C 2 -distance from a single sphere. The corresponding stability 
inequality is obtained with a sharp decay rate. 


1. Introduction 


The aim of this paper is to address a classical question in Differential Geometry, namely the 
characterization of compact embedded constant mean curvature surfaces as spheres - Alexan¬ 
drov’s theorem [4] - in the case of surfaces with constant nonlocal mean curvature. The seminal 
papers [10, 12] have drawn an increasing attention to the geometry of nonlocal minimal surfaces, 
i.e., boundaries of sets 1 l C M n which are stationary for the s-perimeter functional 

p s (n) = f [ . dxdy , 9 , i2 c = ir\i2, 

y ' JaJcic\x-y\ n+2s ' x 5 

corresponding to some value of s € (0,1/2) (see for instance [18, 2, 25, 6, 22, 20, 21, 19]). If 12 
is an open set with smooth boundary and A C M n is an open set, then the condition 

5P s (n)[X] = ^p s ($ t (fi)) = 0, viecn^K"), 

(where denotes the flux defined by the vector-field X ) is equivalent to require the vanishing 
of the nonlocal mean curvature H^(p) of 12 at every point p £ An <912. More in general, we say 
that H ]p : <912 flA-lK is the nonlocal mean curvature of <912 inside A if 

4^($t(12))= / H^{x)X{x)-u x dK~ 1 VleCn^l"). 

at Jan 

Here v x is the exterior unit normal to 12 at x G <912, and "H n_1 denotes the {n — l)-Hausdorff 
measure. 

Whenever <912 is sufficiently smooth (say <912 € C 1,a for some a > 2s), one can show that 
the nonlocal mean curvature of <912 at a point p € <912 is given by 

h ?(p) = ——[ | Xn(x) = xn c (x) ~ Xn(x), (1-1) 

W„-2 J ]R" \X-p\ n+Zs 


where xe denotes the characteristic function of a set E, u n -2 is the measure of the (n — 2)- 
dimensional sphere, and the integral is defined in the principal value sense (see for instance [21, 
Theorem 6.1 and Proposition 6.3]). It is useful to keep in mind that, by means of the divergence 
theorem, the nonlocal mean curvature can also be computed as a boundary integral, that is 


H?(p) = 


1 


L 


(x — p) 


V, 


\n+2s 


x dU™- 1 


( 1 . 2 ) 


SUJn -2 Jdn \x -p\ r 

Our main interest here is describing the shape of open sets 12 having constant, or almost- 
constant, nonlocal mean curvature. In this direction we obtain three main results. 


l 
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The first one is a nonlocal version of the classical Alexandrov’s theorem [4]: 

Theorem 1.1. If Q is a bounded open set of class C 1,2s and H p is constant on dfl, then dVL is 
a sphere. 


In our second result we prove that if Rp, instead of being constant, has just a small Lipschitz 
constant 

\H?(p)-H?(q)\ 


<5 S (II) = sup 

p,q£dn,p^q 


(1.3) 


\p-q\ 

and if dll is of class C 2,a for some a > 2 s, then dfl is close to a sphere, with a sharp estimate 
in terms of To state our result we introduce the following uniform distance from being a 

ball: 

p(Q) = inf { di f im ^ : P e n > B s(p) 


Theorem 1.2. IfQ is a bounded open set with C 2,a -boundary for some a > 2s, then there exists 
a dimensional constant C(n ) such that 


p(Q) < C(n)r} s (£l ), 


(1.4) 


where 


„ zn\ _ diam(fi) 2n+2s+1 x ^ 

— |^| 2 0s(“) ■ 


Moreover there exists rj(n) > 0 such that if r] s (Q) < 77 ( 71 ) then, up 
bi-Lipschitz map F : (0) —>• <90 satisfying 


1 -C(n) y/rj, 
for some dimensional constant C(n ) > 0. 


\x~y\ < \F(x)-F{y)\ < ^1+C(n) \ 


(1.5) 

to rescaling H, we can find a 
-y\ Vi,j/ G dBi(0) (1.6) 


Remark 1.3. Note that both p(H) and r/ s (f)) are scaling invariant quantities. Also, the estimate 
(1.4) is optimal in terms of the exponent of t) s (Q), as it can be easily seen by considering a 
sequence of ellipsoids converging to the unit ball. 

Remark 1.4. If Q is an open set with (7 2 -boundary then (1 — 2 s) Rp —>• H n on dQ as s A 
(1/2) - , where H n is the classical mean curvature of dQ (see [1]). Therefore, because of the 
scaling factor (1 — 2s) one cannot obtain any information from Theorem 1.2 in the limit s —» 
(1/2) - . This is not a drawback of our result, as its local analog is false. Indeed, one can 
construct examples of connected boundaries whose classical mean curvature is arbitrarily close 
to a constant in C 1 topology, but these sets are close (in the Hausdorff distance) to a union of 
tangent spheres of equal radii [8]. 


Both results above are obtained by the moving planes method. Note that the use of this 
method in obtaining stability estimates is well-established in the local case, see for example 
[3, 16, 17] in relation to Serrin’s overdetermined problem and [18] concerning Alexandrov’s the¬ 
orem. Also, this method has already been successfully used in some nonlocal settings to obtain 
symmetry results (see for instance [24, 7] and the references therein). 


Once Theorem 1.2 is proved, we can exploit the regularity theory for nonlocal equations in 
order to obtain a sharp stability estimate in stronger norms. Indeed, by a careful analysis we 
can conclude that dfl is close in C 2 to a sphere with a linear control in terms of rj s (Q), exactly 
as in (1.4). In particular the following result improves the estimate in (1.6), although its proof 
relies on more delicate tools (and actually (1.6) is needed in the proof of this result). 
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Theorem 1.5. Assume that Q is a bounded open set with C 2,a -boundary for some a > 2s, and 
suppose that D has been translated and rescaled so that 

Bi_ Mfl )(0)cfic5i(0). (1.7) 

There exists rj(n,s ) > 0 such that the following holds: If r] s (Q) < r?(n, s) then there is a map 
F : dB\(0) —> M” of class C 2 ' r for any r < 2s, such that F(dB\(0)) = dll and 

||F-Id||cr2,r (SBl( 0)) < C(n,s,r)r] s (n). 

In particular, if r] s (1l) is sufficiently small then H is a convex domain. 


We conclude this introduction by emphasizing that boundaries with constant or almost 
constant mean curvature behave differently in the nonlocal and in the local case, the former 
setting being much more rigid than the latter. Indeed, as proven in Theorem 1.1, even without 
any connectedness assumption a boundary with constant nonlocal mean curvature is a single 
sphere, whereas of course any disjoint union of spheres with equal radii has constant mean 
curvature in the classical sense. Actually, even working only with connected boundaries, a 
significant difference arises at the level of stability. Indeed, as already mentioned in Remark 
1.4, a connected boundary with almost-constant mean curvature may be close to a compound 
of nearby spheres of equal radii (unless one imposes some strong geometric constraints on the 
considered set, like a uniform ball condition [18] or an upper volume density bound [15]). In 
contrast with this picture, as shown in Theorems 1.2 and 1.5 above, uniformly bounded sets with 
almost-constant nonlocal mean curvature must be close to a single ball without the need to any 
uniform control in their geometry. This points out an interesting feature of the nonlocal case, 
namely, the nonlocality of the underlying perimeter functional prevents bubbling phenomena 
(in the limit 5 S ( f2) -A 0). 

We also note that, as it will be apparent from our arguments, Theorems 1.1 and 1.2 hold 
(with different constants and possibly without scale invariant statements) if in the definition of 

one replaces the kernel \z\~ n ~ 2s with k(\z\), where 

A:(t) > 0, t n+2s k{t) + t n+2s+1 \k'{t)\ < C , sup k' < c(t) < 0 , Vt>0. 

(o ,t) 

For the validity of Theorem 1.5, one needs to impose the additional constraint that k[t ) behaves 
as a smooth perturbation of t _ ( n + 2s ) as t —>• 0 + . 

This paper is organized as follows. In section 2 we prove a technical fact about approximat¬ 
ing the nonlocal mean curvature in C 1 with nonlocal “curvatures” coming from smooth kernels. 
Then in section 3 we prove the nonlocal version of Alexandrov’s theorem, while in section 4 we 
address the stability analysis. 

After the writing of this paper was completed we learned that, at the very same time and 
independently of us, X. Cabre, M. Fall, J. Sola-Morales, and T. Weth have proved a result 
analogous to our Theorem 1.1 [9]. 

Acknowledgment: This work has been done while GC and MN where visiting the University 
of Texas at Austin, under the support of NSF-DMS FRG Grant 1361122. In addition, GC 
is supported by a Oden Fellowship at ICES, the GNAMPA of the Istituto Nazionale di Alta 
Matematica (INdAM) and the FIRB project 2013 “Geometrical and Qualitative aspects of 
PDE”. AF is supported by NSF Grant DMS-1262411, and FM is supported by NSF-DMS 
Grant 1265910. 
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2. A TECHNICAL LEMMA 


In order to perform our computations, and in particular to avoid integrability issues, it will 
be useful to work with smooth kernels. We thus consider the approximation K e (x) = <^ e (|x|) 
of the kernel K(x) = — 1 ^ | x |-( n + 2ji ) corresponding to a choice of <p e G ([0, oo)) such that 
tp £ > 0, (p' e < 0, and 


t n+2s v E (t)+t n + 2s+1 \v' £ (t)\<C(n,s) : 
n + 2s 1 


|t 


0J n - 2 t n+2s+1 


as e 


0 + , 


Vf > 0. 


( 2 . 1 ) 


Note that this implies that, as e —>• 0, 


<Pe(t) t 


1 

Un-2 t n+2s 


Vf > 0, 


( 2 . 2 ) 


and both ip £ and ip' £ converge to their limits uniformly on [to, oo) for every fixed to > 0. 
Let us define 

tQ 


Hi 


,{p)= xn(x) (p e {\x - p\) dx, p€dn. 
J R n 


Then, arguing as in [21, Proposition 6.3] we find that 

ll-f^e _ Ilc 0 (an) = 0 > 


(2.3) 

(2.4) 


provided ft is a bounded open set with (7 1 ’“-boundary for some a > 2s. We now prove the 
following technical fact. 


Lemma 2.1. Assume that ft is a bounded open set with C 2,a -boundary for some a > 2s. Then 
H% € C^dft) and H^ £ -A iff in C^dft) as £ ^ 0. 

Proof. Since we already know that H^ e converge to Hg in C° (see (2.4)), it is enough to prove 
that H^ e is a Cauchy sequence in C 1 , that is 

lim \\VHll - VH?J CO(dn) = 0. (2.5) 

To this end we first notice that, by setting 

^e(t) = —p J Te(r) T n_1 dr Vt> 0 

we have 

div (xip e ( |x|)) = ni/} e (\x |) + |x| if' £ {\x\) = <p e {\x\) Vx G M n , 
hence H^ £ can be rewritten as 

H%{p) = ~ 2 I A{\x~p\) {x-p)-v x dH 7 ^ 1 Vped ft. (2.6) 

Jan 

Note that i/y is smooth, it satisfies 

t n+2s ^ e {t) + t n+2s+1 \^' £ {t)\+t n+2s+2 < C(n,s) Vt> 0 (2.7) 

(thanks to (2.1)), and both if e and i\)' £ converge uniformly to their limits on [to,oo) for every 
fixed to > 0 as e —>• 0. 

Now, given p G <912 and e G T p (<9ft) n § n_1 a tangent vector, by the smoothness of ip £ one 
finds 

^ h ?Ap) • e = 2 J (fpe{\x -p\)v x -e + [(* - P) ■ ”x] [(* ~ p) ■ e]) dTtJT 1 . (2.8) 
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Up to decomposing M n = R n 1 x M so that x = (x, x n ) denotes the generic point in K n , and up 
to translating p into the origin 0, we define 

D p = {x £ M n_1 : |x| < p} , Cp = D P x(-p,p), (2.9) 

and we see that the smoothness of dQ implies that, up to a rotation, there exist p > 0 and a 
function / € C 2,a (Dp), with /(0) = V/(0) = 0 and ||/||c 2 ’“(D p ) < £, such that 

Cp n dn = { (x, f{x)) : x £ Dp} = (Id x /)(D p ) 

(by compactness of <90, both p and L are independent of the point p £ dQ under consideration). 
Now, if we set B r = (Id x /)(D r ) for r £ (0,/?), then by the uniform convergence of and 
on [r, oo) we find that 


Jon\B r v 


v x • e + 


V£(W) 


(x ■ v x ) {x ■ e)^ dH 


n— 1 
x 


f (M) Vx ■ e + ^ (x ■ Vx) (x • e)) dU " 1 

dVL\B-r \X\ 


ldn\B- 


0 ( 2 . 10 ) 


as e, rj -A 0. 

On the other hand, having in mind (2.8) and (2.10) and noticing that v x = 

v 1+ i v /( £ )i 

x = ( x , f(x)) £ B r , taking into account that e • e n = 0 for every e £ Tq(< 9U) we see that 


for 


; [ {^e{\x |) 

J ISr 


v x -e + 


V4(M) 


n—1 


= 2 


£ v/-e + 


(x ■ v x ) (x ■ e)^ 

Ve (V\ £ \ 2 + / 


_ [/ — V/ ■ xl (x • e) ) dx, (2-11) 

Vlxp + z 2 1 J v 7 

where / = f{x) and V/ = V/(x). To get a good control on the above quantity, we symmetrize it 
with respect to x by performing the change of variable x i-a — x and then add the two expressions 
(the one with the variable x and the one with —x). In this way we see that the integral in (2.11) 
is equal to 


~ A (VWTJW) (v/(s) • § + Vf(-x) • e) dx 
+ Id x? + fi-xY^j V/(-x)-edx 


+ 


/ 

J D, 


V’e (Vl*l 2 + f(x) 2 ) 
VWFTJW 


+ 


' D, 


(Ve [ V \ x \ 2 + /(®) 2 J 

(\/|^| 2 + /(-^) 2 J 

V \/| x | 2 + /( x ) 2 

\/| x | 2 + /(- x ) 2 


x) — /(—x)] (x ■ e) — [V/(x) • x + V/(—x) • x] (x • e)^ dx 

[/(—x) + V/(—x) • x] (x ■ e) dx. 


Hence, since /(0) = V/(0) = 0 and recalling (2.7), we can find a constant C, depending only 
on n, s, L, such that, for |x| < p, 


|V/(x) • e + V/(—x) • e| < C |x| 1+a , |V/(-x)| < C|x|, |/(x)| < C| 


x 


|/(x) - fi-x )I < C |x| 2 +“, |V/(x) • x + V/(—x) • x| < C'|x| 2+a , 


i>e (v 7 !®! 2 + /(®) 2 ) - V’e (a/I ®| 2 + /(-£) 2 ) 


< c 


1 f(x) 2 - fj-x ) 2 


< 


c 


| x | n + 2 


s—2—ct ’ 
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V4 (Vl^l 2 + f(x) 2 J 

V4 [V\ £ \ 2 + f(-x) 2 ) 

/ r /(t) 2 - /M) 2 

x/\x\ 2 + /Or) 2 

x/\x\ 2 + f{-x) 2 

|£|n+2s+4 


< 


c 


| X | 


thus 


[ (^e{\x\) 
J B r 


v x ■ e + 


04(M) 


(x ■ v x ) (x ■ e)^j dT~L 


71—1 

X 


< C r a ~ 2s , 


( 2 . 12 ) 


where C depends only on n, s and L. Therefore, combining (2.10) and (2.12) we obtain 


limsup|V H^ e (0) • e - Vff s ^(0) ■ e| <Cr 


a—2 s 


£,r )—>-0 


for every r £ (0, p) and any unit tangent vector e € T 0 (dQ) n§" 1 . Hence, by letting r —>• 0 + 
we conclude the proof. □ 


3. Symmetry and the Nonlocal Alexandrov Theorem 


We start by introducing the notation used in exploiting the moving planes method. Given 
e 6 S n_1 , A C R n , and p 6 R, we set 


7T m = {x £ R n : x ■ e = p] 
£/j, = {x £ R n : x ■ e > p} 
Afj = n n 

x’ [l = x — 2 (x ■ e — n)e 
A ',i = {x'v ■ x € A} 


a hyperplane orthogonal to e, 

the half-space on the “positive” side (with respect to e) of 7 t m , 

the “positive” cap of A, 

the reflection of x with respect to n 

the reflection of A with respect to n 

(3.1) 


Now, if is an open bounded (not necessarily connected) set in M n with C 1 -boundary and 
A = sup{x • e : x £ fl}, then for every /i < A sufficiently close to A the reflection with respect to 
7 of the positive cap is contained in H, so it makes sense to define 


A = inf{/7 £ R : (Q^)^ C H for all p, £ (/j,, A)} . 


(3.2) 


In the sequel, given a direction e £ cH?i(0), tt\ and will be referred to as the critical hyperplane 
and the critical cap respectively, and for the sake of simplicity we will set 


x' = x' x = x — 2 (x ■ e — A) e, Vt' = Vl' x = {x r : x £ H} . 


With this notation at hand, we recall from [4] that for every direction e at least one of the 
following two conditions always holds: 


Case 1: dQ' x is tangent to at some point p' £ d£l, which is the reflection in tt\ of a point 
p € dQ x \ tta; 


Case 2: n\ is orthogonal to at some point q £ <9f2 n tt\. 


Both our main results will be based on the analysis of these two possibilities, under the assump¬ 
tion that d s (H) = 0 or that S s (Q) is small, respectively. 


We now prove the following result showing that S s (fi) controls the ^-distance between H 
and fl' (recall that, given two sets E and F , EAF denotes the symmetric difference of the two 
sets, that is EAF = (E \ F) U (F \ E )). Actually, to be able to obtain a sharp stability estimate 
in Theorem 1.2, it will be important to prove a stronger bound on |T2AT2'| when the set H is 
already comparable to a ball of radius 1 (see statement (b) below). 

Proposition 3.1. Assume Q is a bounded open set with C 2,a -boundary for some a > 2s, fix 
e £ § n_1 , and let Q' denote the reflection of Q with respect to the critical hyperplane 7 t\. 
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(a) The bound 
holds with 


|12A12'| < Ci diam(0) n+s+ ( 1 / 2 ) 


Cl = 2 


2 W n —2 


(3.3) 


(3.4) 


n + 2s 

(b) Assume in addition that dist(0,7r,\) < 1/8 and B r ( 0) C 12 C Br( 0) for some radii 
satisfying 

-<r<R<2, R — r > 16 <5 S (12). (3-5) 


Then there exists a dimensional constant C(n) such that 


|i2Ai2'| < c(n) 


Proof. We first prove that 


/ dist(x', n\) dx < — n 2 diam(12) 
Jn aq.' n + 2s 


n+2s+2 


6 s (n). 


(3.6) 


(3.7) 


Without loss of generality we let e = e\. Let us first assume to be in case 1, that is, there exists 
p G <912 a \ tt\ such that p € <912 H <912'. Then 

#?(p) - i/lV) = //.f(p) - tf?'(p) 




1 


-2 \ \x- P \ n + 2s 


dx — 


1 


/ 

1 


1 


| x _p|n+2s 


dx 


(3.8) 


dx , 


w n - 2 7n'\n VI*- 2>l n+2s I*' - p| n+2s , 

where all the integrals are intended in the principal value sense. Since x' = (2A — x\, X 2 , ...,x n ), 

111 r/ \x' — p\ \ n + 2s 


\ X - p \n+2s \ x /_ p \n+2s | x / _ p \n+2s [ \ | x _ p \ 


- 1 


1 


1 + 


4(xi-A)(pi-A )^ 
|x — p| 2 


- 1 


(3.9) 


|x' _p|n+2s 

by the convexity of the function f(t) = (1 + £)( n + 2s )/ 2 — 1 W e get that if x G 12' then 

1 1 2(n + 2s)(xi — A)(pi — A) 2{n + 2s)(x\ — A)(pi — A) 

\x — p\ n+2s \x'— p\ n+2s ~ \x' — p\ n+2s \x — p\ 2 — diam(12) n+2s + 2 

where we used the fact that, by construction, p' G <912 and therefore \x — p\ = \x' — p'\ < diam(12) 
for every x G 12'. Since x\ — A > 0 inside 12' \ 12 and |p — p'\ = 2(pi — A), combining (3.8) and 
(3.9) we find 

2 (n + 2s) 


W ii) > > 


2(pi - A) diam(12) n+2s+2 a; ri _2 7o/\a 
(n + 2s) f 


(xi — A) dx 


n-|-2s+2 


diam(12) 

which proves (3.7) in the first case. 


Wn-2 Jn'An 


\x\ — A| dx . 


We now assume that ir\ is orthogonal to <912 at some point q G <912 n ir\. Thanks to Lemma 
2.1 and (2.3), setting u e (x) = <p e (|x — q |) we have 


vh: 


l (q) ■ e\ = lim VH^Jq) ■ e\ = — lim / Xti(x) Vu e (x) • e\ dx = —2 lim / Vu £ (x) • ei dx 

’ ^oJ r e-+o Jn 
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where we used that f Rn Vu £ = 0. Since X7u e (x) ■ e± = <p' E (\x — g|) is odd with respect 

to the hyperplane {27 = A} (notice that A = qi) and A is the critical value for e\, we find that 
Jantv Vu c • e i = °. hence 

VH^(q) - ei = 2 lim f \ip' E (\x - q\)\ ^ ^ ^ dx. 

e_>0 Jn\si’ F ~ T 

We now observe that 11 \ H' is contained inside the half-space {27 < A} where the function 
non -P os itive, so by ( 2 . 1 ) and monotone convergence we obtain 


VH?(q) • ei = 


2 (n + 2 s) f (x-q)-e i 

w„-2 Jn\n> \x-q \^+ 2 X - 


Since | VH^(q) ■ e±\ < 5 s (fl) and 


(x- q) ■ ei \xi - gi\ 

\x — g | n + 2 «+ 2 — diam(n ) n+2s+2 

we finally get 


\%i - A| 

diam(n ) n + 2s+2 


on 11 \ O' C {x\ < 0} , 


<5 fl (fi) > 


2 (n + 2 s) 

diam( 01 ) n + 2 s+ 2 a ; n _2 



A| dx = 


(■n + 2 s) 


diam(n) n + 2 s + 2 w n _ 2 Jo'AO 


L 


\x\ — A| dx , 


which completes the proof of (3.7). 


We now prove (a). For this it is enough to combine (3.7) with Chebyshev’s inequality to get 

|{x G flAfi' : dist(x, 7 r^) > 7 )! <- n 2 diam(n) n+ 2 s+ 2 (5g(n), 

7 n + 2 s 

that together with the trivial bound 

|{x € flAfi' : dist(x, Tt\) < 7 }| < 27 diam(fl) n_1 , 

gives us (3.3) choosing 7 = 2 (n+ 2 s) diam(fl) s+ ( 3 / 2 ^ y/5 s (Al). 

If we know in addition that dist( 0 , 7 TA) < 1/8 and that B r { 0 ) C 11 C Br{ 0 ) for some radii 
satisfying (3.5), then we can use the stronger bound 

|{x€l2Afl': dist(x, 7 Ta) < 7 }| < C{n) 7 (R — r) V 7 <l/ 4 , 

so (3.6) follows by choosing 7 = \J. □ 

We now deduce Theorem 1.1 from Proposition 3.1. 

Proof of Theorem 1.1. We begin by noticing that, thanks to the regularity theory developed 
in [ 6 ] (see in particular the proof of [ 6 , Theorem 1]), C l,2s domains with constant nonlocal 
mean curvature are actually C°°, so Proposition 3.1 applies. In particular, since by assumption 
<5 S (11) = 0, Proposition 3.1 implies that 11 is symmetric in any direction. 

Since the barycenter b of H belongs to every axis of symmetry and every rotation can be 
written as a composition of reflections, we have that H is invariant under rotations, which implies 
that <911 is a collection of concentric spheres centered at b. To show that till is just one sphere, 
we apply again the method of moving planes in an arbitrary direction: if 911 is not connected 
then the critical hyperplane must be a hyperplane of symmetry and cannot contain b, which is 
a contradiction. Hence 911 must have a single connected component, i.e., 911 is a sphere. □ 
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4. Stability 


Before proving Theorems 1.2 and 1.5 we first show the following lemma stating that if <5 S (12) 
is small then, up to a translation, all critical planes from the moving planes method pass close 
to the origin. Again, as in Proposition 3.1, it will be important to show a stronger bound when 
11 is comparable to a ball of radius 1. 


Lemma 4.1. Let 12 be an open bounded set of class C 2,a for some a > 2s with 


diam(n) n+s+ ( 1 / 2 ) 


v^M< min {^}y 


n + 2s 
8co n -2 


(4.1) 


and suppose that the critical planes with respect to the coordinate directions ir ei coincide with 
{xi = 0} for every i = 1, ...,n. Also, given e € § n_1 , denote by X e the critical value associated 
to e as in (3.2). 

(a) The bound 

|A e | < Csv 7 ^) (4.2) 

holds with 

diam(12) n+ ' s+ ( 3 / 2 ) 

C 2 — 4 (n + 3)-j^j- C \, 

where C\ is as in (3.4). 

(b) Assume in addition that dist(0,7r A ) < 1/8 and B r ( 0) C 11 C Br{ 0) for some radii 
satisfying (3.5). Then 

|A e | < C%n)y/6&i)VR=T (4.3) 

for some dimensional constant C*(n). 


Proof. We first prove (a). To this aim, we define H° = {—x : x € H} and set 

C* = C x diam(H) n+s+(1/2) , (4.4) 

where C\ is defined as in (3.4). Then, since H° can be obtained from H by symmetrizing it with 
respect to the hyperplanes {xi = 0} = 7T ei for i = 1,..., n, applying Proposition 3.1 with respect 
to the coordinate directions we obtain 

|HAll 0 1 <nCiy/5 a (n). (4.5) 

Now, to prove (4.2) we assume that X e > 0 (the case X e < 0 being similar). We first note that 

A e = sup{x ■ e : x € H} < diam(H). (4-6) 

Indeed, if A e > diam(H), then x ■ e > 0 for every x € H, and thus |HA!2°| = 2|11|, which 
contradicts (4.5) and (4.1). This said, we denote by fl 7 the reflection of H about the critical 
hyperplane ir\ e , and deduce from Proposition 3.1 that 

|UAH'| < Ciy/6 s (n). (4.7) 

Now, recalling the notation 12^ = 11 n = H n {x ■ e > p,}, it follows by (4.7) (which tells us 
that H is almost symmetric with respect to n\ e ) that 

|^A e | > - C? y/SjJl). (4.8) 

Since H is almost symmetric about 0 by (4.5), using the notation = {—x : x € £\ e } we see 
that (4.8) gives 

| 12 n£° I = |12°n£ Ae | > |H a | — |HAH°| > H - (n + 1)^^, 

2 
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which together with (4.8) implies 

\{x £ fl : — A e < x ■ e < A e }| < (n + 2)C*\J <5 s (n). (4-9) 

In other words, by combining the almost-symmetry of II with respect to 0 and to 7r Ae we have 
shown that II has small volume in the strip {|x • e| < A e }. Since {A e < x ■ e < 3A e } is mapped 
into {\x ■ e\ < A e } by the reflection with respect to 7r Ae , exploiting again (4.7) and (4.9) we get 

|{x £ Q : A e < x ■ e < 3A e }| = |{x £ fi' : lx • el < A e }| 

__ (4.10) 

< |{s € n : |x • e| < A e }| + |DAD'| < (n + 3)C^5 S (Q). 


Define now 

rrik : = |{x £ ft : (2 k — l)A e < x ■ e < (2k + 1)A e }| , k > 1 , 

and notice that, by the moving planes procedure, the set D n (seen as a subset of M™^ 1 ) is 
included inside D n tt^ whenever A e < /r' < /i. In particular the function [i i —> 'H Tl ~ 1 (Q, n i r^) is 
decreasing on (A e , A e ), hence m^ is a decreasing sequence and (4.10) gives us 

mk < mi < (n + 3)Cj \J d s (D) V k > 1. 

Recalling that D C {x • e < A e }, combining this last estimate with (4.9) and letting ko be the 
smallest natural number such that (2fco + l)A e > A e we get 


I^A e 


ko 

n {A e < x ■ e < A e }| < ^2 m k A 

k=\ 




hence (thanks to (4.6)) 

|^A e | A e < (n + 3) C* diam(fl) y/SJSTj. 
Since |OA e | > |D|/4 (by (4.8) and (4.1)), recalling (4.4) we get (4.2). 


To prove (b) it suffices to observe that, under the assumption that dist(0,7TA) < 1/8 and 
-B r (0) C D C -Br(O) with r,R satisfying (3.5), we can repeat the very same proof done above 
but using (3.6) in place of (3.3) to obtain (4.3). □ 

We now prove Theorems 1.2 and 1.5. 

Proof of Theorem 1.2. Step 1: proof of (1.4). Up to a translation, we can assume that the 
critical planes with respect to the coordinate directions TT ei coincide with (x* = 0} for every 
i = l,..., n. 

Notice that, since p(Vt) < 1 and i /rjjJV) = dldlll( ^j -one can directly assume 

that (4.1) holds. Moreover, setting 

r = min |x| , R = max |x|, (4.11) 

x£dfl xddSl 

it is enough to control R — r (as it gives an upper bound on p(Q)). 

Let x,y £ dQ be such that |x| = r and \y\ = R. Assuming without loss of generality that 
x / y, we consider the unit vector 

y-x 
I y-x\' 

and let tt\ b denote the corresponding critical hyperplane. We notice that y is closer than x to 
the critical hyperplane n\ e , i.e., 

dist(x,7r A J > dist(y,7r A J . (4.12) 

Indeed, since x = y — te with t = |x — y\, the method of moving planes implies that the critical 
position can be reached at most when y' (the reflection of y with respect to 7r Ae ) is tangent to 
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x, which corresponds to the equality case in (4.12), while in all the other cases strict inequality 
holds. Thus, by (4.12) and the fact that e is parallel to y — x we get 

R - r = \y\ - \x\ < 2 dist( 0 , 7 r A J = 2 |A e | (4.13) 

that combined with (4.2) implies that 

R-r < 2CZy/5 s (Q) = 16 (n + 3) t/ 2 ^ n ~ -^(O). (4.14) 

V n + 2s 

Now, since all the quantities involved are scaling invariant, we rescale 12 so that R = 1, and we 
assume without loss of generality that 

R-r> 16(5,(17) 

as otherwise (1.4) trivially holds. In this way it follows from (4.14) that (3.5) holds provided 
77 ,( 12 ) is small enough. Also, thanks to (4.2) we see that dist(0, 7 r Ae ) < 1/8 for all e € § n_1 if 
(5,(12) (or equivalently 77 ,( 12 )) is sufficiently small. 

Hence, this allows us to combine (4.13) with (4.3) to get 

R-r < 2C*(n)yj5 s (p) VR^r, 

that is 

R — r < 4 C* (n) (5, (12), (4.15) 

which proves (1.4). 

Step 2: a quantitative Lipschitz bound on (912. We want to show that if 77 ,( 12 ) < 77 ( 71 ) for some 
dimensional constant 77 ( 77 ,), then <912 is Lipschitz-flat with a uniform bound. 

Since all the quantities involved are scaling invariant, we assume as at the end of step 1 that 
R = 1 so that 

B r ( 0) C 12 C .Bi(O) 

with 

1 — r < C(n) rj s (12) (4.16) 

(by (4.15)), and then prove (1.6) for 77 ,( 12 ) small enough. 

To this end, it is enough to show that there exists a dimensional constant M = M(n ) such 
that, for any x € (912 and y 6 dBi_ Mris ^( 0) such that the “open” segment (x,y) is contained 
outside Bi_ Mr)a ( q)( 0), then (x,y) C 12. Indeed, this means that for any x € <912 we can find a 
uniform cone of opening 7 r — C^/rjjJV) with tip at x and axis parallel to ^ which is contained 
inside 12, and this implies that <912 is locally the graph of a Lipschitz function satisfying (1.6). 

Now, to prove the latter fact, assume by contradiction that there exist x € <912 and y € 
9Bi_Mr] a (n)( 0) for which there exists a point z € (x,y) H 12 c . Set e = ZZy\ an< 4 notice that, 
since z € 12 c , it follows that the moving planes method has to stop before reaching z, that is 
A e > z- e. Now, since (x,y) C Hi(0) \ H 1 _ M7?s(Q) ( 0) and y € <9-Bi_ M 7 ?s(Q )( 0), we have y ■ e > 0. 
Hence, since z — y is parallel to e and z € 12 c C B r ( 0) c we get 

A e > (z - y) ■ e + y ■ e > (z - y) ■ e = \z - y\ > M 77 ,( 12 ) - (1 - r). 

On the other hand (4.2) gives 

0 ( 77 ) 77 ,( 12 ) > |A e | 

(recall that 1 < diam(12) < 2 and s € (0,1)), which leads to a contradiction to (4.16) provided 
M = M{n ) is large enough. □ 
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Proof of Theorem 1.5. Our goal here it to exploit the results from Theorem 1.2 to get closeness 
to a ball in higher norms. For this, we need to show that our assumptions on iLp imply that 
dPt is smooth with some quantitative bounds depending only on rj s (Q). Hence, we first formulate 
the following regularity criterion that is implicitly contained in [14] (recall that definition of C r 
and D r from (2.9)). 

Given n > 2, s,£ G (0,1/2), and (3 G (0,2s), there exist positive constants e = e(n,s,l,/3) 
and C* = C*(n,s,£, /3) with the following property: Let E be an open set with C 2 -boundary 
such that for some L > 0 it holds 

\\Hf\\c° { dE)<Li \ Br ^ n n E \ € V y G dE ,r < £. (4.17) 

Ld n 1" 

If 0 G dE and r < t are such that 

B r (0) n dE C {x € M n : \x n \ < er} , Lr<e , (4.18) 

then there exists u € C' 1 ’^(D r / 2 ) such that 

C r / 2 n dE = (Id x u)( D r/2 ), 

with 

l|Vu||co(D r/2 ) + 7-/3 [V«]c°.?(D r/2 ) < C* (- r ( r) + Lr^j . 


Step 1: uniform C 2,7 bounds on dPt. We show that the regularity criterion stated above applies 
with E = fi. 

Since, by the definition of y(fl), the radii 1 — 2p(fX) and 1 must be optimal for the inclusion 
(1.7) to hold, we can find points p\ € D ^^ 1 _ 2 p(n)(0) and p 2 G f2n3Hi(0). Hence, it follows 
by the inclusions (1.7) and (1.1) that 

< Hs 1 ^, H?( P2 ) > 

and because the Lipschitz constant of is bounded by S s ( fi) < C rj s (Q) and 

\Hs h - 2p(n) - Hf l \ < Cp(n) < CpsiQ) 


(by (1.4)), we deduce that 

\\H?-H*\\ Looidn) <CVs(to). (4.19) 

Notice now that the uniform Lipschitz estimate provided by Theorem 1.2 implies that the 
density estimates in (4.17) hold. Thus, provided y s (H) is small enough, (4.17) holds with L = 
2 H^ 1 and for some £ = £(n) > 0. At the same time we can find r = r(n) > 0, depending on 
dBi (0) only, such that if x G dB\ (0) then 


B 2r (x) n 5Hi(0) C |y G M n : (y - x) 


< 


er 


}> 


Lr < e. 


(4.20) 


Hence, assuming that y s (H) is small enough in terms of r, by (1.4) and (1.7) we can ensure that 

z 

\z I 


B r {z) PI dQ C |y G M n : (y — z) ■ — < erj VzG d£l , 


(4.21) 


and applying the regularity criterion stated before we obtain that, for any z G d£l, there exists 
a uniform neighborhood such that, in a suitable system of coordinates, dQ is given by the graph 
of a function u z : D r —» R with 


ll“allch/J(D r/2 ) < C(n,s,(3). 

Now, choosing j3 arbitrarily close to 2s and exploiting the fact that G C°’ 7 (clH) for every 
7 G (0,1) together with the higher regularity theory by [6, Section 3], we obtain that 
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for any r < 2s. 

Step 2: dQ is C 2 -close to a sphere linearly in r] s (£l). By the previous step we know that there 
exists a map / : dB i(0) -A M of class C 2,T for any t < 2s satisfying 

ll/llc 2 ^(dBi(o)) < C(n,s,r) 

and such that dfl = {y + f{y) y : y 6 <9i?i(0)}. Notice that, by (1.7), 

ll/II.L oo (aBi( 0 )) < C(n) rj s (Q) , (4.22) 

so we deduce by interpolation that for any C, <2s there exists an exponent a(C) > 0 such that 

ll/llc 2 .c(aB 1 (o)) < C(n, s, C) rjs(n) a ^ . (4.23) 

This implies in particular that dQ is C 2 -close to a sphere, so is convex for rj s (Q) sufficiently 
small. We now want to show that (4.23) is still valid if we replace a(C) with 1, which will prove 
the theorem with F(y) = y + f(y) y. 

For this, we write the nonlocal mean curvature in terms of / starting from (1.2): in this 
way, since any point x E can be written as y + f(y)y with y E <9i?i(0), by the area formula 
we get that, at the point p = q + f(q ) q E dQ, 

flfM = — / i l + t <T v - q - t /[ t \U ■ (»- ^rrn) (> + m ” 1 <™r‘- 

su n - 2 j dBl (o) \y + f(y)y-q-f(q)q\ n+2s v 1 + fiy) J y 

To simplify the notation we define the vector-field v q (y) := y + f(y) y — q — f(q) q, so that the 
above expression becomes 

v q {y) 


h?(p) = -2— I 

SUl n -2 Jd 


-2 J dBl ( 0 ) \v q (y)\ n+2s 


[y{l + f{y)) n (1 + f{y)Y 1 -(! + /(?))" 1 ) 




Now, noticing that the normal to <9-Bi(0) at y is equal to y itself, by the tangential divergence 
theorem (see for instance [23, Theorem 11.8]) we get (notice that the classical mean curvature 
of dBi(0) is n — 1) 

1 f .. / v q (y) 


Hs(p) = 


+ 


s(n- l)cu n -2 JdB 1 ( 0 ) 


dfw 


K(y)l 


1 


L 


v q (y) • y 


n-\-2s 


(1 + f(q)) n l dU n y -\ 


(i + /(y)) B_1 -(i + /(3)) n_1 


dH 


n— 1 


SUJ n -2 J 9Bl ( 0) \Vq(y)\ n+2s 
Since di vt{v) = n — 1 and S/rf(y) • y = 0 we have 

v g (y) 


\Vq(y) I 


dlw 
So, computing 

and denoting by n y : R n 

v q (y) ■ V T \v q (y)\ = 


n-\-2s 


(n — l)(l + /(y)) c^ v g(y) • VrK(y)| 


Ml/) I 


n-\-2s 


— (n + 2 s)- 


My)l n+2s+1 


V TVqiy) = (l + f(y)) Vry + v r /(y) ® y 

-A M n the orthogonal projection onto y we get 
Vq{y) • V TVqiy ) ~ Vq{y) 

\v q {y)\ 

(! + f(y)) hyVq(y )| 2 - (l + /(<?)) ((<? - y) ■ V T /(y)) (My) • y) 


My) I 


Thanks to the elementary identity 


(y - q) ■ y = 1 - q ■ y = -j\v - q \ 2 


( 4 . 24 ) 
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we see that 


and 


\ 7r y v g(y)\ 2 = (! + f(.q)) 2 W y q\ 2 = (1 + f{q)) 2 { 1 +yq) ^ ^ , 


wg(y) • y = f(y) - /(g) + (i + /(g)) 


|y - q ? 


Hence, setting for simplicity 


r f{y,q) = {l + f{y)) n 1 - (l + /(<?))" 1 


and combining all these formulas, we obtain 




\p) = —— I 

SUJ n -2 Ja 


i + f(y) 


-2 Zasbo) k(y)l n+2s 


r /(y,g)dW 


71 — 1 
2/ 


n + 2s 


(i + f(y)) (i + y ■ q) (i + /(g)) 2 1 y - q I s 


2s (n - l) w n _ 2 /asi(o) M^)l n+2s+2 


r/(y,g)d^ 


n—1 


+ 


+ 


n + 2s 


s (n - 1) ca „_2 


f (! + /(g)) ((g ~ y) • Vr/(z/)) (/(g) - /(g)) ^ ^ mn - X 

Jai ^■' iIti+2s+2 ’ 1 


/aei(o) 


n + 2s 


L 


\v q (y)\ n+2s+2 

(i + /(g)) 2 ((g - y) ■ v r /(y)) |y - g| 2 


2s UJ n —2 
Noticing that 


2s (n - 1) w „_2 /asi(o) 

+ -^(!+ /(,))"-'/ 
•S w n _2 Jq 

(! + /(?))"/ 
Jd 


? l -|- 2 s -|-2 


\Vq{y) I 

/(2/) ~ /(g) ~,n-l 

ob i( o) K(y)l n+2s y 


Ff(y,q)dK 


n— 1 

y 


\y - q? 


dU 


n —1 


sb i( o) k(y)l n+2s y 


l + yg = 2--|y-g| 2 , 


the above expression can be rewritten as 




*(p) = — / 

SUJn-2 Jd 


1 + /(z/) 


r/(y, g) dJ-L 


n —1 
2/ 


-2 ./asi(o) |w g (y)| n+2s 

n + 2s f (i + /(y)) (i + /(g)) 2 |y - g| 2 


s (n — 1) ca n _2 /asi(o) 


L 


My) I 


?i-|-2s-|-2 


r/(y, g) dU 


n —1 


+ 


+ 


n + 2s 


s (n — 1) w n _2 /asi(o) 


L 


n + 2s 


L 


(i + /(g)) ((g - y) • Vt /(z/)) (/(g) - /(g)) 

My)l n+2s+2 

(i + /(g)) 2 ((g - z/) • v T /(z/)) |y - g| s 


r /(y,g)d^ 


n—1 

y 


2s (n - 1) W n _2 /aBi(O) |u g (?/)| rl + 2s + 2 


r/(j/,g)dw: 


n—1 

y 


j-Tjn—l 

SB,(0) l 1, e(y)l n+2s 5 


+ —(1+ /W)- 1 / ^ /(,) 

«W„-2 JdB 1 (0) \Vq(y)\ 


+ 


+ 


2s W„-2 
1 


4s LO n —2 /aRi(O) 


L 


SBi(O) Mz/)l' 

(i + /(z/)) (i + /(g)) 2 lz/-g| 4 

|t>q(y)| n+2s+2 


r/(y,g)rfH 


71 — 1 

y 


(4.25) 
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We now notice that, since 

r f(y,q) = (n- l)[/(y) - /(g)] (1 + P(f(y), 

with P(t, s) a polynomial of degree n — 2 which vanishes at t = s = 0, the first five terms in the 
right hand side above can be written as 


I {f(y)-f(q))K(y,q)dn n y - 1 
JdBi (0) 


IdB i(0) 

where the kernel K(y,q ) behaves like a C 1,T perturbation of the -fractional Laplacian on 
R n_1 : more precisely 


K{,J ’«) = • 


(4.26) 


where Gf : dB\(0) x cXBi(O) — > M is a C' 1,r -function (depending on /) which satisfies 
l|G , /llc 1 ^(asi(o)xasi(o)) < C' ll/llc 2 ' T (aBi(o)) Vt€ [0,2s). 

We now subtract the value of the above expression in the right hand side of (4.25) at / = 0 
(which corresponds to the case of the unit sphere) to get 

H?(F(q)) - = - [ {f(y)-f(q))K(y,q)d'H n y - 1 +g(q), (4.27) 

JdB, for 


where F(q) = q + f(q) q and 
1 


g(q) = 


+ 


2 S U) n —2 


(1+/(?))’ 


/aBi(o) 


\y-q\ 2 

\v q {y)\ n+ ^ s 


dU™- L - 




l 


aBi(o) 


— r/77 

| y _ ^|ra+2s—2 U7 W 


n—1 


-/ 

G n-2 Jai 


4s W n _ 


(i + /(y)) (i + f{q)) 2 \y - q\ 4 

\n+2s+2 


r f(y,q)dn 


71—1 

y 


is a C 1 function satisfying 

H5llL°°(aBi(0)) < C II/IIl°°(9Bi(0))) llfflie 1 (aBi(O)) < C ll/llcraBiCO))- 

Since K is a C l perturbation of the 1 ^ 2g -fractional Laplacian, applying [11, Theorem 61] locally 
in charts (using a cut-off function) we deduce that 

ll/llc 1 >' r (9Bi(o)) < C( n , s,T)(\\f\\i J °o(QB 1 (o))+\\g\\L°°( d B 1 (o))P\\H^oF—HB 1 \\ L oo( 0 Bl ( O ^ Vr < 2s. 

Also, differentiating (4.27) we can apply the same result to the first derivatives of / (see for 
instance [6, Section 2.4] for more details on how this differentiation argument works) to get 

H/||c 2 . t (9Bi(0)) < C(n, S,T)^||/|| C i( 9S i( 0 )) + ||ff||c' 1 (3Bi(0)) + ll^°-P 1 -^f 1 ||c 1 (9Bi(0))) Vr < 2s - 
Notice now that by (4.19), the definition of <5 s (fi), and the fact that ll^Ucirnsj^)) < C, we have 

\\H?oF-H*\\ cHaBim <C6 9 (n). 

Hence combining all these estimates and recalling (4.22), we conclude that 

< C(n, s,t) rj s (Q,) Vr<2s, 


□ 


c 2 ’ r (aBi(o)) < C(n, s,r)(5 s (H) + 


as desired. 
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